NUHEP-TH-94-20 

Aug 1994 



S-WAVE AND P-WAVE B c MESON PRODUCTION AT HADRON 
COLLIDERS BY HEAVY QUARK FRAGMENTATION 

KINGMAN CHEUNG* 
Department of Physics, Northwestern University, Evantson, Illinois 60208, U.S.A. 

. . . and 

TZU CHIANG YUAN 
Davis Institute for High Energy Physics, Department of Physics, 
o I University of California, Davis CA 95616, U.S.A. 

<D 

in 

Os ■ ABSTRACT 



in 



We compute model-independently the production rates and transverse mo- 
mentum spectra for the B c mesons in various spin-orbital states (n 1 So, n 3 Si, 
n 1 Pi, and n 3 Pj (J = 0, 1, 2) ) at hadron colliders via the direct fragmentation 
of the bottom antiquark and via the Altarelli-Parisi-induced gluon fragmen- 
tation. Since all the radially and orbitally excited states below the BD flavor 
threshold will decay, either electromagnetically, hadronically, or a combination 
of both, into the pseudoscalar ground state l 1 ^, they all contribute signifi- 
cantly to the inclusive B c meson production. 

Oh! 

Ph! The next and the last family of B mesons to be observed will be the B c (bc) made 

up of one charm quark and one bottom antiquark. Like the J/ip and T quarkonia, 
dynamical properties of B c can be predicted reliably by using perturbative QCD, in 
contrast to the heavy-light mesons. In the limit m c /m& — > 0, the B c system enables 
us to test the heavy quark symmetry and to understand the next-to-leading terms in 
the heavy quark effective theory in the applications to the heavy-light B mesons. In 
addition, the production rates for different spin-orbital states also help us to understand 
the spin symmetry breaking effects. Phenomenologically, B c mesons can be used to 
analyze the mixing of the B® — B® without ambiguity by tagging the charge of the 
lepton in the decay B+ -> 5° + £ + u or B~ + £~v 

Calculations on the production of B c mesons at e + e~ colliders were previously 
performed. But the calculation for the production at hadronic colliders is rather te- 
dious until Braaten and Yuantl pointed out that the heavy quarkonium production at 
the large transverse momentum region is dominated by heavy quark fragmentation. 
The fragmentation of a heavy quark into a heavy-heavy-quark bound state essentially 
involves the creation of a heavy quark-antiquark pair, which tells us that the process 
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should be hard enough to be calculable in perturbative QCD (PQCD). Explicit cal- 
culation of the production rates and the transverse momentum spectra was performed 
in Ref.,0 which was based on the^ direct b antiquark fragmentation functions D b ^ Bc (z) 
and D b ^ B , (z) calculated in ReiB It was founda that the Altarelli-Parisi-induced gluon 
fragmentation g — > B c also contribute significantly to the iotal production. Since the 
P-wave fragmentation functions have just been completed p it is natural to include all 
the S-wave and P-wave contributions to calculate the inclusive production rate. 

The direct b — > B c fragmentation function and the induced g — > B c fragmentation 
function are coupled by the following Altarelli-Parisi equations: 

d dy dy 

^-^--Ds-Bc = / —P b ^ b {z/y,fi) D- b ^ Ba (y,fi) + / —P- b ^ g (z/y,ti) D g _ Bc (y^) , 
c/i Jz y Jz y 



(1) 

l*a~ D s->B c (z, H) = — p g ^b(z/y,^)D- b ^ Bc (y,fi)+ —P g ^ g (z/y, fj) D g ^ Bc (y, //) . 

Oil Jz V J z V 



d f 1 dy r 1 dy 

(2) 

where Pi_>j can be approximated by the usual massless Altarelli-Parisi splitting func- 
tions. Similar equations can be written down for the 3 Si, Pi, and 3 Pj (J = 0,1,2) 
states. The boundary conditions for the coupled equations are D g ^ Bc (z, //) = for 
fi < 2(m b + m c ) and D b ^ Bc (z, fiQ = m b + 2m c ), which is the heavy quark fragmentation 
function calculated to the order of a 2 s at the initial scale /iq by PQCD. Expressions for 
the initial fragmentation functions in .S'-wave states can be found in RefB and those 
for P-wave states can be found in Ref.u 

Numerically integrating the coupled equations with the above boundary con- 
ditions, we obtain the direct b antiquark fragmentation functions and the induced 
gluon fragmentation functions for the S'-wave and P-wave states at any arbitrary scale 
A 4 > ^o- The inputs to these fragmentation functions are the quark masses m b and m c 
and the nonperturbative parameters associated with the wavefunction of the bound 
state. These nonperturbative parameters can be calculated within the framework of 
the potential modelsJ3 For the two S'-wave states there is only one nonperturbative 
parameter, which is the radial wavefunction P(0) at the origin. The P-wave fragmen- 
tation functions have two nonperturbative parameters associated with the color-singlet 
and color-octet mechanisms. Two of the P-wave states ( x Pi and 3 Pi) mix to form two 
physical states because they have the same quantum numbers. The two physical states 
are denoted by |1+) and The mixing and further details can be found in Ref.B 

The calculation of B c meson production is simplified by factorizing the whole 
process into a short- distance process of producing the heavy quark and a long-distance 
process, which is the fragmentation of the heavy quark into the B c meson. The differ- 
ential cross-section for the B c meson in the x So state is given by 

da(B c (p T )) = J2 J h/ P {xi,v)fj/v{ x 2,v) [dcr(ij -> b(p T /z)X, fi) D b ^ Bc (z, n) 

+ da(ij -> g(p T /z)X, D g ^ Bc (z, //)] , (3) 



where i,j denote all the possible initial partons. The first term in the square bracket 
is the direct b fragmentation contribution and the second term is the induced gluon 



Fig. 1. The differential cross sections da / dpx{Bc) versus the transverse momentum px(B c ) 
of the B c meson for different spin-orbital states at the Tevatron. 

fragmentation contribution. Similar expressions can be written down for the 3 Si and 
the P-wave states. In the above equation, the factorization scale /i is chosen to be of 
the order of Pt/z to avoid large logarithms in da, while the large logarithms in D(z) 
can be summed up by evolving the D(z) according to the coupled equations in Eqs. ([[]) 
and (0). 

In our calculation the factorization scale \i in Eq. (fj) is chosen to be 

= vH, 9 + m b ' ( 4 ) 

where px h is the transverse momentum of the fragmentating parton. p? h is related 
to pt(B c ) by PT lg = Pt(B c )/z. This scale is also used for the parton distributions 
and the running coupling constant a s . Explicitly we used m h = 4.9 GeV, m c = 1.5 
GeV, and CTEQ20 for the parton distributions. The a s (Q) is evaluated at 1-loop by 
evolving from the experimental value a s {m z ) = 0.118 by a s (Q) = a s (m z )/(l + ((33 — 
2nf)/67r)a s (mz) log((5/mx)),_where n/ is the number of active flavors at the scale Q. 
We included gg — > bb, gb — > gb, and qq —>■ bb as the hard subprocesses for the inclusive 
b production, and gg — > gg, gq(q) — > gq(q), and qq — > gg for the inclusive g production. 
The fragmentation functions at the scale /i are obtained by solving Eqs. (|l|) and (0) 
with the boundary conditions mentioned above. 

The transverse momentum spectra for different spin-orbital states are shown in 



Table 1. The integrated cross sections in pb for the B c mesons in various spin-orbital states, 
with the acceptance cuts in Eq. (|5|), at the Tevatron. |1+) and are the two physical 

P-wave states resulted from the mixing of the 1 P± and 3 Pi states. 
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Fig. |TJ, in which the contributions from the direct b fragmentation and the induced 
gluon fragmentation have been added, with the acceptance cuts 

p T (B c ) > 10 GeV and \y(B c )\ < 1 (5) 

on the B c mesons. The curves in Fig. |TJ are for n = 1 states. The radially excited n = 2 
states can be calculated similarly with the corresponding nonperturbative parameters. 
We present the integrated cross sections for n — 1 and n = 2 states in Table [l|. 

Since the annihilation channel for the decay of the excited B c meson states is 
highly suppressed relative to the electromagnetic and hadronic transitions, all the ex- 
cited states below the BD threshold will decay into the ground states by emitting 
photons or pions. Thus, they all contribute to the inclusive production. Adding all the 
contributions shown in Table |I|, we have a total cross section of 1.2 nb, which implies 
about 1.2 x 10 5 B r mesons for 100 pb 1 at the Tevatron. This number should almost 
represent the total inclusive rate, except for a small contribution from the D-w&ve 
states. 

Thus, we have presented the so far most complete B c meson production via heavy 
quark fragmentation, including 5-wave and P-wave contributions, at the Tevatron. 
The B c meson can be detected via the decays into J/if> + X, in which the J /if) is fully 
reconstructed by the leptonic decay. If X is a charged lepton, then the event has a very 
distinct signature of three charged leptons coming out from the same displaced vertex. 
If X can be fully reconstructed, together with the reconstructed J/ if) the B c meson can 
be fully reconstructed. 
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